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Abstract 

We discuss the computation of the grand canonical partition sum describing hot mat- 
ter in systems with the Higgs mechanism in the presence of non-zero conserved global 
charges. We formulate a set of simple rules for that computation in the high-temperature 
approximation in the limit of small chemical potentials. As an illustration of the use of 
these rules, we calculate the leading term in the free energy of the standard model as a 
function of baryon number B. We show that this quantity depends continuously on the 
Higgs expectation value 0, with a crossover at ~ T where Debye screening overtakes 
the Higgs mechanism — the Higgs vacuum "melts" . A number of confusions that exist in 
the literature regarding the B dependence of the free energy is clarified. 
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Screening of gauge fields by the Higgs mechanism is one of the central ideas both 
in modern particle physics and in condensed matter physics. It is often referred to as 
"spontaneous breaking" of a gauge symmetry, although, in the accurate sense of the 
word, gauge symmetries cannot be broken in the same way as global ones can. One might 
rather say that gauge charges are "hidden" by the Higgs mechanism. For example, in the 
standard model at zero temperature, particles are characterized by values of the electric 
charge, but not the isospin and the weak hypercharge separately.^] 

At high temperatures, another mechanism for screening becomes operative — the Debye 
screening of charges in plasma. The Debye screening is not associated with "breaking" 
or hiding of gauge charges; on the contrary, it is due to the motion of conserved charges. 
Thus, in the standard model, while all elementary excitations (particles) at zero tem- 
perature (T = 0) have only one conserved gauge charge with respect to electroweak 
interactions, it could be different for excitations at T 7^ 0. 

In general, we can envisage a competition between the two screening mechanisms. The 
mass of vector bosons generated by the Higgs mechanism is of order gcj), where g is the 
gauge coupling and is the temperature dependent Higgs expectation value. The electric 
mass due to the Debye screening is of order gT. So, when T <^ (p, we expect the vacuum 
classification of particles to be more or less intact; when T ^ 0, the system forgets that 
the gauge symmetry is "broken" . 

At T ~ a crossover occurs, which can be called a melting of the Higgs vacuum. This 
regime may or may not be associated with a phase transition. A gauge theory may have 
no gauge-invariant order parameter, and in that case a phase transition does not have to 
occur [U, 3. In fact, it was recently shown that the electroweak phase transition changes 
into a smooth crossover when the mass of the Higgs boson becomes sufficiently large 0. 
In this sense, the minimal standard model is similar to the liquid-gas system: there is no 
true distinction between the phases, and the line of first-order phase transitions ends in 
a critical point. 

The melting of vacuum, associated with a crossover from one screening mechanism to 
the other, admits a fully gauge-invariant characterization using the renormalized average 
of (p^(p, where is the corresponding Higgs field. In the extreme high-temperature limit, 
in a weakly coupled theory, 

{(P^(P)r = i^— + {coTiections), (1) 

where p is the number of components of if, and the corrections are due to interactions. As 
the temperature lowers, the leading term starts to deviate from (at zero temper- 

ature it is Temperatures at which the deviation becomes significant (of the order 

of the leading term itself) mark the melting crossover. 

The above considerations show that the couplings of elementary excitations to gauge 
fields (the charges) at high temperature can be different from those at T = 0. This effect 

^ Although we will use the standard model as an illustration, our results are not restricted to that 
case and can apply, for instance, to grand unified theories. 
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can enter calculations of thermodynamic quantities in systems with the Higgs mecha- 
nism, which contain, in addition, densities of some conserved global charges. Therefore, 
a treatment of such systems requires some care. As we will see below, the use of kinetic 
equilibrium requirements in terms of zero-temperature particle excitations, sometimes in- 
voked 0, |, ^ in these circumstances, may, and in fact does, lead to wrong results at 

Although the only application we discuss here will be for the electroweak sector, it is 
of some interest to sketch a general formalism that would allow one to deal with any case 
when both the Higgs and the Debye screening mechanisms are present. 

In general, a system will have a number of conserved global charges, which we denote 
by Ni. These are gauge- invariant quantities that allow for the introduction of chemical 
potentials, /Xj, in the usual way. We thus can define a gauge- invariant grand partition 
sum as a Euclidean functional integral: 



exp{-Q/T)= Jv<l>exp ~ J dTd^x{LE f^^^^ 



(2) 



where $ denotes generically all of the fields of our system; the integral over r is restricted 
to the finite interval from to /5 = 1/T, with the usual periodic ( ant iper iodic) conditions 
for bosons (fermions). The Lagrangian Le includes gauge-fixing terms and ghosts. The 
thermodynamic potential is a gauge-invariant function of fii and T. 

Special attention should be paid to integrals over the Euclidean temporal components 
of the gauge fields. The integration over their zero-momentum modes, the set of which 
will be denoted by A4, enforces the conditions of neutrality of our system with respect 
to all gauge charges. These conditions are simply a consequence of the Gauss constraint, 
integrated over the whole three-dimensional (3d) space.0 

It is convenient to imagine that the functional integration in (H) is performed in two 
steps. In the first step, we integrate out all modes except for A4, and the Higgs fields ip. 
Those are integrated out at the second step. The reason for this two-step procedure is 
that ^4 and ip can develop expectation values, which we have to take into account. Notice 
that this general method of calculation of Cl makes no use of the notion of a particle or 
an elementary excitation in plasma. 

The result of the first step is an effective Euclidean action Q 



exp(-f2/r) = yWexp - J dTd^x{LE -J2 l^i^i) 



(3) 



where the prime on the integral shows that the ^4 and ip integrations are omitted, is a 
function of /Xj, T, Ai and (p. One may notice that iAi play the role of chemical potentials 
for the corresponding gauge charges. In general, f2 is a gauge non-invariant quantity; the 
integration over A4 and ip will convert it into the gauge-invariant potential Q. 

Let us now restrict ourselves to weakly-interacting cases when we can compute Q by 
the loop expansion. We will also assume that the expectations values of A4 and ip can be 

^For definiteness we consider the 3d space as a compact manifold with large volume which we eventually 
send to infinity. 
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Figure 1: Tadpole diagram producing terms linear in chemical potential. 



chosen x- and r- independent. The calculation can be done in any desired order; here we 
will need only tree and one-loop terms. At the tree level, the effective action is the sum 
of the covariant derivative term for and its tree potential, 



where in the covariant derivative we set d^ip = 0; V is the total 3d volume. One-loop 
terms can be classified according to powers of A4. Terms independent of A4 add to Vq 
to produce the usual temperature-dependent potential for the scalar field. Then, there 
are linear terms (tadpoles, see Fig. 1), quadratic terms (Debye masses) as well as higher 
powers. We will assume that the chemical potentials for the global charges are small, 
/ij/T -C 1. In that case, the expectation values of A4, are proportional to fii, and the 
terms beyond the quadratic ones are suppressed. In addition, the tadpoles and the Debye 
masses can be expanded in /ij/T. For the tadpoles, the expansion starts with terms linear 
in /ij, and for the Debye masses we can, in the leading order, neglect /ij altogether. 

A particularly simple set of rules arises in the important case when, in addition to the 
above limits, the masses of all particles are much smaller than T. Here the results can 
be simply formulated using the spectrum of the high-temperature, "unbroken" phase. In 
that case, the one-loop effective action is obtained as follows. To each particle species of 
the high-temperature limit we associate a chemical potential 



where n[ are the global charges of that species and gj its gauge charges; / labels the 
species. Note that only A4,, corresponding to the mutually commuting generators of 
the gauge group, should be taken into account. These generators can be simultaneously 
diagonalized, and their eigenvalues define what we mean by the charges in eq. 

■^For the hypercharge, our convention is gy — lli/'^- 



no/v = \DM^ + Vo{v). 



(4) 



(5) 
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The leading terms in the one-loop effective action are 



Q,/v = [Vt{^) - v,{^)] -t'y: (6) 

where rji = 1/6 (77/ = 1/12) per single helicity bosonic (fermionic) particle-antiparticle 
pair. The complete effective action needed for our purposes is the sum of the tree action 
(|) and the one- loop action d). 

As an illustration of the use of these rules, we present here a calculation of the leading 
term in the free energy of the standard model (SM) as a function of the baryon number. 
It generalizes a similar calculation done in ref. 0, where the only fermionic degrees of 
freedom were the leptons of the first generation. The B-dependent part of free energy 
plays an important role in the theory of anomalous electroweak B non- conservation. In 
our previous calculation of this quantity [0], we pointed out that at sufficiently high 
temperatures, one may have to use the spectrum of the "unbroken" phase even below 
the phase transition. Here we extend our previous result and present the leading term 
in F{B) for all temperatures T ^ ■m^r{T). In particular, we demonstrate a crossover in 
F{B) at T ~ 0. 

As is well known, there are Uf strictly conserved global charges in the SM with massless 
neutrinos, namely 

N, = —B- Li, (7) 

where i = l,...,nf is the index of fermionic generation. In addition to them, there are 
two gauge charges associated with the underlying SU(2)xU(l) symmetry, which commute 
with the Hamiltonian, the global charges defined by eq. (|^), and each other. These are 
the hypercharge Y, 



i=l 



1 - 4- 2- 



and the third component of the weak isospin T3, 

= e,,,4Qo + l^f^r'^Do- Do r^)f + ^ E [Q^lor^Qi + ^i7orUJ . (9) 

i=i 

Here Qi and Li are the left quark and lepton doublets, respectively; Ui, Di are the right 
quark fields, and Ei are the right leptons. 

According to our general rules, we have to introduce nj chemical potentials /ij to con- 
served numbers (^, and compute the effective action ^(^44, B^, fii, 4>), where (p^ = (0, -^). 
To make the equations somewhat more transparent we assume here that the asymmetries 
in charges iVj are degenerate in flavour. This assumption is harmless for the leading-order 
calculation. It would be inadequate if we wanted to include corrections in leptonic masses, 
which are important in some cosmological scenarios |]^. This assumption sets all /Zj equal 
to each other and thus replaces them with a single chemical potential for B — L. For 
future purposes, though, we will keep separate chemical potentials for B and L. The 



4 



^')/T^ (12) 



condition of equilibrium with respect to the anomalous B-non-conservation, jJ,B + I^l = 0, 
can be imposed later. 

The leading one-loop terms of the high-temperature expansion of the effective action, 
for small chemical potentials, are0 

n{Al B,, /X,, 0) = ^m2(T)02 + iA0^+ 
l<P'{gAl - g'B,f + \mlAlAl + ]^m'lB,B, (10) 

where 

< = + (11) 

2 fins 

m'(r) ^ 4,,4 + r^(iA + ^,^ + i-,«+i„.). (13) 

and (7t is the Yukawa coupling of t-quark. Terms proportional to /i^,yu| are the standard 
contributions of the fermionic chemical potentials to fi, while the terms linear in chemical 
potential come from the tadpole diagrams of Fig. 1. 

The equation (|l^) forms a basis for determining the equilibrium properties of the hot 
and dense electroweak plasma at small chemical potentials. As a first example, let us find 
the equilibrium value Bq of the baryonic number at fixed value of the strictly conserved 
charge B — L = X . That equilibrium value is obtained from the solution of the system 
of equations: 

(neutrality of the system with respect to gauge charges); 

— = -7^ = ^ (15) 

(definitions of average baryonic and leptonic charges); the equilibrium condition 

1^3 + 1^1 = 0, (16) 

and the normalization condition B — L = X . We get 

^ 4(2n/ + n,)mf, + 8(2 + 2n/ + n,)m^^(0) 

° (22n/ + 13n,)m|) + 2(26 + 24n/ + 13n,)m|i,(0) ^ ' ' ^ ' 



We omit 0^ term in the effective potential which is not essential for the present discussion. 
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Figure 2: Dependence of the equilibrium value of the baryon number (in units of B — L) 
on the expectation value of the Higgs field (in units of temperature) for Uf = 3, Ug = 1. 



where mw[(f)) = \g(j). Even if there is a phase transition between the symmetric and 
Higgs phases, in the Higgs phase with </> -C T (the Debye mass is much larger than the 
Higgs-induced vector mass), we still get the same result as in the symmetric phase R 



In the opposite limit ^ T we get 



4(2ny + 



22n/ + 13n, 



-{B-L). 



Bo 



4(2 + 2n/ + n,) 
26 + 24n/ + 13n, 



[B-L) 



(19) 



which coincides with the result of refs. @, |^, |, § • In the intermediate region, the value 
of the baryon number interpolates between the two limiting cases, see Fig. 2. 

For cosmological applications, the relation between Bq and X should be taken at the 
moment of the freeze-out of the sphaleron processes. For a strongly first-order phase 
transition, when after the transition 0/T >1.2-1.5, the freeze-out coincides with the 
transition itself, and to a good accuracy one can use the result for the symmetric phase, 
eq. (18). For a weakly first-order transition, a second-order transition, or a crossover, one 
should use the full eq. ( |17|) with (p/T ~ 1.2-1.5. (That corresponds to m^(0)/m|) ~ 0.2- 
0.3.) Looking at Fig. 2, we see that for these values of 0/T, the equilibrium value of B 
is approximately in the middle between its two limiting values. 

Although the numerical difference between the limiting values is small, we believe 
that this calculation is worth-while because it elucidates the physics responsible for the 
crossover between the high- and low- temperature regimes. What matters here is the 
relation between cj) and T, and not T and Tc, as stated in refs. |Q, |], or T and mw{(t>) 
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Figure 3: Dependence of the coefficient k on the expectation value of the Higgs field (in 
units of temperature) for n/ = 3, = 1. 



as stated in ref. 0. The error of the equilibrium reaction analysis of refs. |^ is in 

the use of the zero-temperature spectrum at high temperatures, which led the authors of 
these references to conclude that the chemical potential of the Higgs boson is equal to 
zero. In our notation, that would correspond to g'B/^ — gA\ = 0. As can be seen directly 
from eq. (0) with the effective action ([T0|), this condition indeed holds at T = 0, but not 
at any T 7^ 0. 

Another quantity of interest is the free energy of the system at fixed baryonic number 
B, not equal to Bq, which is used in the computation of the baryon erasure rate due to 
sphalerons [0. It is defined as 



F{B) = n + fiBB + fiLL, 



(20) 



where baryonic and leptonic chemical potentials are to be found from eqs. (|T^,|T5|) but 



the equilibrium condition (16) is now not imposed. The result is 



F{B) 



Bo? 



VT2 



const, 



where 



3[(22n/ + 13n,)m% + 2(26 + 24n/ + 13ns)ml,{(^)] 
4n/[(5ra/ + 3n,)m|) + (12 + lluf + 6ns)m^(0)] 



For small 0/T we have 



3(22n/ + ISris 
Anfi^Uf + 3n5 



(21) 



(22) 



(23) 
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which for = 1 coincides with eq. (4.23) of [0. For large 0/T: 

^ ^ 3(26 + 24n; + 13n,) 

2raj(12 + llnj + 6n,) ^ ^ 

The dependence of k on 0/T is shown in Fig. 3. 

In conclusion, we have presented a general way of dealing with conserved global quan- 
tum numbers in finite-temperature gauge theories with the Higgs mechanism. We have 
formulated a simple set of rules for computing the grand canonical partition sum in the 
high-temperature limit. Using these techniques, we have computed the leading term of 
the high-temperature expansion of the free energy of the standard model, as a function 
of baryon number. We have shown that this quantity depends continuously on the Higgs 
expectation value 0. The crossover in the free energy is associated with a "melting" of 
the Higgs vacuum at ~ T, when the Debye screening overtakes the Higgs mechanism. 
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